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I 

REQUIRED  THRESHOLD  SETTINGS  AND  SIGNAL-TO-NOISE 
RATIOS  FOR  COMBINED  NORMALIZATION  AND  OR-ING 

INTRODUCTION 

When  weak  signals  of  unknown  strength  and  location  have  to  be 
detected  in  the  presence  of  noise  of  unknown  and  varying  level, 
it  is  necessary  to  make  estimates  of  the  intensity  of  the 
interfering  background.  These  estimated  (noise)  levels  are  then 
compared  with  that  for  a  candidate  signal  level  and  location  for 
purposes  of  making  statements  about  the  likelihood  of  signal 
presence  or  absence  in  that  particular  data  segment  under 
investigation.  Here,  we  will  investigate  the  performance 
capability  of  such  a  normalizer,  in  terms  of  the  false  alarm  and 
detection  probabilities,  and  determine  the  thresholds  and  input 
signal-to-noise  ratios  leguired  to  attain  these  probabilities. 

Furthermore,  when  large  amounts  of  multichannel  data  have  to 
be  processed  with  limited  computational  facilities  in  reasonable 
or  short  intervals  of  time,  it  is  necessary  to  resort  co 
shortcuts  or  data  reduction  in  order  to  avoid  overload-  One 
possible  approach  is  to  employ  or-ing,  in  which  only  the  largest 
of  a  set  of  quantities  is  retained  for  further  data  processing 
and  decision  making. 

Finally,  in  practice,  it  is  often  necessary  to  utilize  both 
normalization  and  or-ing  together.  This  combination  of  nonlinear 
processors  requires  a  resetting  of  the  thresholds  that  would  have 
been  appropriate  for  use  of  the  normalizer  alone.  Here,  we  shj=*‘'i 
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investigate  all  three  situations,  namely  normalization,  or-ing, 
and  a  combination  of  normalization  and  or-ing,  in  terms  of  the 
probabilities  and  thresholds  stated  above. 

When  the  size  of  the  nonlinear  transformation,  whether  it  be 
normalization,  or-ing,  or  both,  is  changed,  the  required 
thresholds  will  have  to  be  changed  if  the  previously  realized 
(false  alarm)  probabilities  are  to  be  maintained.  For  example, 
suppose  we  had  been  or-ing  8  random  variables  and  decided  to 
change  the  or-ing  size  to  accept  16  random  variables  instead,  for 
purposes  of  further  data  reduction.  Then,  the  required  threshold 
settings  would  have  to  be  modified  to  maintain  specified  false 
alarm  probabilities,  as  would  the  required  input  signal-to-noise 
ratios  for  specified  detection  probabilities.  This  maintenance 
of  probabilities  under  a  change  of  size  of  transformation  will  be 
investigated  here. 

At  the  .^aire  time  that  the  size  of  a  nonlinear  transformation 
is  changed,  it  may  be  desired  to  subject  its  output  to  a 
different  number  of  thresholds  than  ’-'ere  utilized  previously. 

This  possibility  is  allowed  and  analyzed  in  addition. 

The  particular  physical  situation  considered  here  is  that  of 
multiple  simultaneous  beamformer  outputs,  each  with  large  banks 
of  narrowband  filters  subject  to  envelope-squared  detection.  For 
Gaussian  noise  inputs,  the  probability  density  function  of  each 
of  these  device  outputs  is  exponential.  Furthermore,  when  a 
Gaussian  signal  is  also  present  at  the  input  of  one  of  these 
narrowband  filters,  the  corresponding  probability  density 
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function  of  the  detected  output  is  still  exponential,  but  with  a 
level  governed  by  the  signal-to-noise  ratio  at  that  particular 
filter  output.  This  scenario  will  be  the  mainstay  of  the 
analysis  here. 

The  physical  motivation  for  this  study  is  to  be  able  to  set 
requantization  levels  on  a  display,  in  order  to  achieve  constant 
marking  density  independent  of  the  signal  processing  parameters 
such  as  the  amount  of  or-ing  and  normalizer  size.  Such  displays 
occur  in  active  as  well  as  passive  sonar  systems. 
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DEFINITIONS  OF  FUNCTIONS 

The  random  variables  (RV)  for  noise-only,  for  0  <  m  <  M, 

are  independent  and  identically  distributed  (IID)  with  common 
probability  density  function  (PDF)  P^i^)  and  cumulative 
distribution  function  (CDF)  C^(X),  where 

X 

C^(X)  =  Prob(x^  <  X)  =  j  du  p^(u)  .  (1) 

—  00 

The  corresponding  characteristic  function  (CF)  of  RV  is 

=  I  exp(i^u)  Px(n)  ,  (2) 

where  j.ntegrals  without  limits  are  over  The  moments  and 

cumulants  of  order  j  of  RV  x^  are  denoted  by  and 

respectively.  The  inverse  function  to  CDF  C^(X)  in  (1)  is 
denoted  as  C^(P);  that  is, 

if  P  =  C^(X),  then  X  =  C^(P)  for  0  <  P  <  1  .  (3) 

As  an  example,  consider  narrowband  filter  outputs  for  which 
p^ ( u )  =  exp(-u)  for  u  >  0  , 

C^(X)  =  1  -  exp(-X)  for  X  >  0  , 

f^(l)  =  (1  -  iJ;)"^  for  all  I  , 

C^(P)  =  -  ln(l  -  P)  for  0  <  P  <  1  , 

Pjj(j)  =  j!  ,  =  (j-l)J  for  j  ^  1  .  (4) 
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The  scaling  of  RV  in  (4)  has  been  taken  such  that  its  mean 
is  1.  This  is  done  solely  for  notational  convenience;  it  will 
not  affect  the  probabilities  realized  herein  nor  the  required 
signal-to-noise  ratios.  However,  it  does  influence  the  threshold 
settings  calculated,  which  would  have  to  be  scaled  for  a 
different  input  noise  level. 
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THRESHOLD  RESETTING 

A  collection  of  IID  RVs,  for  0  1  m  1  M,  is  subject  to  a 

nonlinear  transformation  f,  yielding  output 

y  =  f(XQ,Xj,...,x„)  .  (5) 

For  example,  a  normalizer  is  characterized  by 


y  = 


+  x 


M 


Mil, 


(6) 


while  an  or-ing  device  yields 

y  =  max(Xj^,X2, . . .  ,Xj^)  ,  M  i  1  .  (7) 

A  combination  of  a  normalizer  and  or-ing  device  will  require  a 
more  general  formulation;  then  we  would  use 


X 


k 


M 


w, 

ko 


M 


for  1  <  k  i  K  , 


y  =  max(Xj^,X2, . . .  ,Xj^)  ,  (8) 

where  we  need  two  parameters,  M  and  K,  and  are  K(M+1)  IID 

RVs.  Of  course,  then  the  K  RVs  are  also  IID. 

Let  Cy(Y)  be  the  CDF  of  RV  y  at  the  output  of  general 
nonlinearity  f  in  (5).  Choose  threshold  Yj  such  that  specified 
cumulative  probability  is  realized  there;  that  is. 
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Pi  =  Cy(Yi)  ,  or  Yi  =  Cy(Pi)  ,  (9) 

where  is  the  inverse  function  to  C^.  1“^!  is  the  false  alarm 

probability  at  threshold  Yi*  Also,  choose  additional  thresholds 
lYi.j  such  that  Yi  <  Y2  <  •••  <  Y^,  for  a  total  of  T  thresholds, 
with  the  largest  one  being 

Yt  =  ^1  Ay  ,  where  dBy  =  10  logiQ(Ay)  (10) 

is  a  specified  decibel  value.  Then  take  the  remaining  thresholds 
according  to  equal  spacing  rule 


Y 


t 


for  1  1  t  i  T  . 


(11) 


The  normalized  thresholds  for  RV  y  are  defined  as 


Xt  ' 


s 

for  1  1  t  1  T  . 


(12) 


where  Py  and  Cy  are  the  mean  and  standard  deviation  of  RV  y. 

Then  compute  the  cumulative  probabilities  realized  at  these 
thresholds  {^1.},  namely 


P  =  C  (Y  )  for  1  i  t  £  T  , 
u  y  ^ 


(13) 


and  print  out  M,  dBy,  Py,  Cy,  T,  fY^j,  {Y^l,  {P^}*  The  false 
alarm  probabilities  are  {1-P^}. 
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An  alternative  choice  is  to  space  the  thresholds  (Y^j  equally 
in  decibels.  That  is,  defining 

=  10  logj^Q(Y^)  for  1  1  t  i  T  (14) 

as  the  threshold  in  decibels,  we  take  as  given  in  terms  of  Yj^ 
and  we  take 

+  dBy  .  (15) 

The  intermediate  decibel  thresholds  are  then  selected  according 
to  the  equal  spacing  rule 


D 


t 


for  1  1  t  1  T  . 


(16) 


The  power  thresholds  can  then  be  evaluated  as 

D  /lO 

Y^  =  10  for  1  i  t  i  T  .  (17) 

The  cumulative  probabilities  at  these  latter  power  thresholds 
follow  from  (13)  as  before. 

Another  possibility  is  to  simply  set  the  thresholds  according 
to 

Yt  =  ^y  +  <^y  Xt  ^  t  1  T  ,  (18) 

where  normalized  thresholds  {Y^)  are  preset  constants  determined 
by  the  user.  In  this  latter  case,  probability  in  (9)  would 
not  be  realized  at  initial  threshold  Yj  =  /Jy  +  Y ^ .  In  any 
event,  the  desired  printouts  are  the  quantities  listed  under 
(13)  . 
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Finally,  we  could  set  the  thresholds  such  that  preset  values 
of  probabilities  {P^l  are  realized  for  all  1  1  t  i  T.  That  is, 
solve  (13)  for  the  thresholds  according  to 

1  i  t  1  T  .  (19) 

This  amounts  to  setting  T  different  false  alarm  probabilities 
{l-Ptl- 
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ALTERNATIVE  SIZE  TRANSFORMATION 

Now  consider  the  new  RV  z  obtained  by  changing  the  parameter 
value  from  M  to  N  in  the  given  nonlinear  transformation  in  (5): 

2  ~  i ( Xq » /  •  •  • » Xjj )  •  (20) 

N  can  be  larger  or  smaller  than  M.  Let  the  CDF  of  RV  z  be  C^(Z). 
We  now  choose  new  thresholds,  {Z^}  for  1  £  t  1  T,  such  that  the 
probabilities  {P^l  in  (13)  are  maintained  for  RV  z;  that  is,  we 
choose  thresholds  (Z^)  for  RV  z  in  (20)  such  that 

Cz(Zt)  =  Pt  '  or  Z^  =  Cjj(Pt)  for  1  i  t  i  T  .  (21) 

These  thresholds  (Z^l,  to  be  employed  for  RV  z,  will  not 
necessarily  have  equal  spacing,  as  did  the  thresholds  in 

(11),  for  example,  for  RV  y.  The  normalized  thresholds  for  RV  z 
are 

^t  ■  ^z 

Z.  =  — - -  for  1  i  t  i  T  ,  (22) 

— t  CT  '  ' 

Z 

where  u  and  a  are  the  mean  and  standard  deviation  of  RV  z. 

'^z  z 

Then  print  out  N,  dB  ,  p  ,  o  ,  T,  {Z.  },  {^+.1,  {P^.},  where 

dB^  =  10  logjQ(Z^/Zj)  .  (23) 
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DIFFERENT  NUMBER  OF  THRESHOLDS 


If  RV  z  is  to  be  subject  to  a  different  number,  U,  of 
thresholds  than  RV  y  was,  it  is  not  always  reasonable  to  try  to 
maintain  the  set  of  T  probabilities  realized  in  (13).  (U 

can  be  larger  or  smaller  than  T. )  One  alternative  is  to  maintain 
the  edge  probabilities  and  according  to  (21),  thereby  de¬ 
termining  values  for  and  Z^.  Then  choose  a  different  complete 
set  of  thresholds  {Z^}  for  RV'z  according  to  equal  spacing  rule 

Z'  =  Z,  +  -n - i-^(u  -  1)  for  1  i  u  i  U  .  (24) 

We  can  then  evaluate  the  cumulative  probabilities  at  these 
latter  threshold  values  as 

P'  =  C  (Z' )  for  1  £  u  i  U  .  (25) 

u  z'  u'  '  ' 

Of  course,  P^  =  P^  and  Py  =  P^,  since  Zj  =  Z^  and  Zy  =  Z^.  This 
also  means  that  dB  is  still  given  by  (23). 

It  must  be  noted  that  this  change  in  philosophy,  namely  to 
maintain  only  edge  probabilities  Pj^  and  P,p,  will  not  reduce  to 
the  earlier  results  when  we  set  U  =  T  here.  The  new  thresholds 
in  z,  given  by  (24),  are  equally  spaced  even  if  U  =  T,  whereas 
the  former  thresholds  given  by  (21)  are  not. 

Print  out  N,  dB^,  o^,  U,  (Z'j,  {Z'),  {P'j,  where  {Z^}  are 

the  normalized  thresholds 


Z' 


u 


-  P, 


a 


z 


for  1  1  u  i  U  . 


(26) 
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SIGNAL-TO-NOISE  RATIO  REQUIREMENTS 

Let  US  return  to  the  original  nonlinear  transformation  (5) 
characterized  by  parameter  M.  The  CDF  of  RV  y,  with  signal 
present  in  just  one  of  the  input  RVs  {Xjjjl/  is  denoted  by  Cy(Y;R), 
where  R  is  the  input  signal-to-noise  ratio  (SNR)  in  that 
particular  RV  containing  a  signal.  If  we  want  this  new  CDF  to 
realize  a  specified  probability  value  Sj^  at  all  the  thresholds 
lY^l  in  (11),  the  required  SNRs  {R^l  must  satisfy 

^1  "  for  1  1  t  1  T  .  {2/, 

From  (13)  and  (27),  we  know  that 

=  Cy(y^)  =  Cy(Y^;0)  for  1  i  t  <  T  .  (28) 

Therefore,  specified  probability  Sj^  in  (27)  must  satisfy 

S^  <  for  1  i  t  i  T  ,  (29) 

in  order  that  nonnegative  SNRs  {R^.1  can  result  as  solutions  to 
(27).  That  is,  each  cumulative  distribution  value  must  be 
decreased  from  to  Sj,  meaning  that  each  exceedance  probability 
has  been  increased  from  false  alarm  probability  value  1-P^  to 
detection  probability  value  Pdj  =  1-Sj^.  The  actual  determination 
of  CDF  Cy(Y;R)  will  have  to  wait  until  after  we  have  specified 
and  analyzed  the  nonlinear  transformations  (5)  of  interest,  for 
the  case  of  noise-only  present. 
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If  several  detection  probabilities,  like  Pd^  =  .5,  Pd2  =  .9, 
Pda  =  .99,  are  specified,  there  will  be  a  set  of  equations  like 
(27)  for  each  case,  namely 

1  -  =  Sj  =  Cy(Y^;R^)  for  1  1  t  i  T  ;  j  =  1,2,3  .  (30) 
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STATISTICS  OF  NORMALIZER 


Let  RV  y 
from  I ID  RVs 


be  obtained  by  means  of  a  normalization  procedure 
0  i  m  1  M,  according  to  transformation 


y 


+  X 


(31) 


where  we  assume  that  x^^  2  0  for  all  m.  The  average  RV,  a,  in  the 
denominator  of  (31)  has,  respectively,  CF  and  PDF 


Pa^^^  n  J  exp(-iu^)  [f^(?;/M)]^ 


(32) 


The  CDF  of  RV  y  in  (31)  is,  using  the  statistical 
independence  of  Xq  and  a, 

Cy(Y)  =  Prob(y  <  Y)  =  Prob(XQ  <  Ya)  = 

00  Yt  <® 

=  J  dt  Pg(t)  J  du  p^(u)  =  J  dt  Pg(t)  C^(Yt)  .  (33) 

0  0  0 

The  corresponding  PDF  of  RV  y  is,  upon  differentiation, 

oe 

Py(u)  =  J  dt  t  p3(t)  Px(ut)  .  (34) 

0 

The  moments  oi  RV  y  can  be  found  from  a  couple  of  forms: 
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Py(j)  =  =  J  du  Py(u)  = 

0 

=  xj  a~^  =  ^^(j)  ^^(-j)  =  J  du  ^  J 

Convergence  of  the  last  integral  in  (35)  may  not  occur  for  larger 
values  of  j;  that  is,  due  to  the  division  in  (31),  RV  y  nay  not 
have  finite  higher-order  moments. 


EXAMPLE 

The  example  presented  in  (4)  is  utilized  here.  From  (4) 
and  (32),  the  CF  and  PDF  of  average  RV  a  in  (31)  are 


fjl)  = 


(1  -  i£;/M) 


M 


(36) 


and 


Pa(u) 


m"  u”-*  exp 1 -mu)  for  u  >  0 

(M-1) !  °  u  ?  u  , 


(37) 


respectively.  The  CDF  of  RV  y  then  follows  from  (33)  and  (4)  as 


Cy(V)  =  j  dt  «" 
0 


=  1 


-M 


for  Y  >  0  . 


(38) 


The  PDF  of  RV  y  is  therefore 
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Py(»)  *  (l  b) 


-M-1 


for  u  >  0  . 


(39) 


The  inverse  function  to  CDF  Cy(Y)  in  (38)  is 


C  (P)  =  M  (1  -  P)  -  l]  for  0  <  P  <  1  . 


(40; 


The  j-th  moment  of  RV  y  is  given  by 


,(j)  =  y^  =  J  du  Py(«)  =  J  du  (l  + 


=  B(i  +  l,M-i)  =  r(j-H)  r(M-j)  ^  j.!,M3. 


where  we  used  (39)  and  (1;  3.194  3  and  8.384  1] 


M 


2M' 


y  ~  M-1  '  y  (M-2 ) (M-1 ) 


'  2 
^  (M-2)(M-1) 


M' 


-1 

mJ 


-M-1 


for  j  <  M  ,  (41) 

In  particular, 

for  M  >  3  .  (42) 


The  condition  M  >  3  is  necessary  for  RV  y  to  possess  a  finite 
variance.  We  now  have  all  the  quantities  required  for 
application  in  the  previous  section  on  threshold  resetting. 
For  M  =  ®,  f^[^)  in  (36)  equals  exp(i£;);  that  is,  RV  a  in 

a 

(31)  is  equal  to  the  constant  1.  This  corresponds  to  no 
normalization  and,  in  fact,  (31)  reduces  to  y  =  x^.  Also, 
(38)  -  (41)  reduces  to  (4),  as  expected. 
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CHANGE  IN  SIZE  OF  NORMALIZER 

If  we  now  change  from  size  M  to  N  in  normalizer  (31),  we 
obtain  RV  z  as  considered  in  (20)  and  sequel.  Its  CDF  follows, 
by  similarity  of  form  to  (38),  as 

C^(Z)  =  1  -  (l  -r  I)  for  Z  >  0  .  (43) 

Its  inverse  function  is 


c^(p: 


=  n[( 


-1/N 


-  1 


Reference  to  (42)  also  reveals  that  the  mean  and  standard 
deviation  of  RV  z  are 


for  N  >  3  . 


(45) 


The  new  thresholds  are  given  by  (21)  and  (22). 
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STATISTICS  OF  OR-ING  DEVICE 


Let  (x  j,  1  1  m  1  M,  be  IID  RVs  with  conunon  PDF  p  (u)  and  CDF 
C  (X).  The  CDF  and  PDF  of  the  maximum  RV 

X  '  ' 


y  =  max{Xj,X2,  .  .  .  ,Xj^}  , 


(46) 


yielded  by  an  or-ing  device,  are  then 


Cy(Y)  =  [C^(Y)]^  , 


Py(u)  =  M  (C^(u)]”  ^  p^(u)  , 


(47) 

(48) 


respectively.  The  inverse  function  to  CDF  C^  in  (47)  is 


~  ~  (  l/M) 

C^(P)  =  C^  P^^” 
y  A  V 


for  0  <  P  <  1 


(49 


Here,  again,  C  is  the  inverse  function  to  CDF  C  of  RV  x  . 

^  X  X  m 

The  CF  of  RV  y  follows  from  (48)  as 


f  (^)  =  I  du  exp(i^u)  M  [C^(u)]^  ^  • 


(50) 


The  moments  of  RV  y  may  be  obtained  from  (48)  in  the  form 


,(j)  =  J  du  u^  M  [C^(u)]”  ^  • 


(51) 


Alternatively,  the  cumulants  can  be  obtained  by  a  power  series 
expansion  of  the  natural  logarithm  of  CF  fy(J;)  in  (50). 
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EXAMPLE 

We  again  use  the  results  given  in  (4).  Substitution  into 
(50)  yields 


^  (1)  =  J  du  exp(i?;u  -  u)  M  [1  -  exp(-u)] 


M-1 


M  I  dt  t  ( 1-t)”  ^  = 


rd-i^;)  r(M+i) 
r(M+i-ij;) 


M 


TTfi  -  ^1 

[  mj 


m= 


-1 


(52) 


where  we  let  t  =  exp(-u)  and  used  [1;  8.380  1  and  8.384  1].  This 
is  a  very  compact  form  and  is  easily  computed  numerically,  if 
necessary.  This  result  illustrates  that  RV  y  has  the  same 
statistics  as  RV  w  defined  by 


M 


w  = 


m=0 


jn 

m 


1*2 


M*M 


(53) 


In  order  to  determine  the  cumulants  of  RV  y,  we  consider 


In  fy(J:) 


M 


m=l 


m=l  k=l 


(54) 
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The  cumulants  follow  immediately  as 


Xy(j) 


M 


for  j  2  1  . 


(55) 


In  particular,  the  mean  and  variance  of  RV  y  are 


=  Xy(l) 


(56) 


respectively.  It  is  seen  that  the  mean  of  RV  y  increases  without 

limit,  in  fact,  logarithmic  with  M.  On  the  other  hand,  the 

2 

variance  saturates  at  n  /6,  meaning  that  the  standard  deviation 
Oy  of  RV  y  cannot  exceed  n/6^  =  1.283. 

The  CDF  of  y  follows  from  (47)  and  (4)  as 

Cy(Y)  =  [1  -  exp(-Y)]”  for  Y  >  0  ,  (57) 

with  corresponding  inverse  function 

Cy(P)  =  -  ln(l  -  for  0  <  P  <  1  .  (58) 

For  M  =  1,  there  is  no  or-ing,  and  (46)  reduces  to  y  =  Xj^. 
Also,  (52),  (55),  (57),  and  (58)  reduce  to  (4),  as  expected. 
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CHANGE  IN  SIZE  OF  OR-ING  DEVICE 

If  we  change  from  size  M  to  N  in  or-ing  device  (46),  we 
obtain  RV  z  considered  in  (20)  and  sequel.  Its  CDF  is,  by 
similarity  in  form  to  (57), 

C^(Z)  =  [1  -  exp(-Z)]^  for  Z  >  0  .  (59) 

The  corresponding  inverse  function  is 

C^(P)  =  -  ln(l  -  for  0  <  P  <  1  .  (60) 

Reference  to  (56)  reveals  that  the  mean  and  variance  of 
RV  z  are  given  by 

N  .  ,  N  T 

Pz  =  ZZ  H  '  <^2  =  m  -i  •  (61) 

^  n=l  ^  ^  n=l  n^ 

The  new  thresholds  for  RV  z  are  given  by  (21)  and  (22). 
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STATISTICS  OF  NORMALIZER  AND  OR-ING  DEVICE 


Here,  we  consider  a  set  of  K(M+1)  IID  RVs  {w.  )  subject  to 


both  normalization  and  or-ing,  according  to 


=  M  ^  '^km  '  = 

m=l 


for  1  i  k  i  K  , 


y  =  max(Xj^,X2/ • . .  ,Xj^) 


The  statistics  of  the  normalization  portion,  namely  for 

1  1  k  i  K,  were  previously  considered  in  (31)  -  (35)  for  the 
general  case,  and  then  specialized  to  an  example  in  (36)  -  (42). 
Also,  the  analysis  of  the  or-ing  portion  was  conducted  in 
(46)  -  (51)  and  then  specialized  to  an  example  in  (52)  -  (58). 

We  will  rely  heavily  on  those  results  in  order  to  minimize  the 
presentation  in  this  section. 


EXAMPLE 


We  presume  that  all  the  input  RVs  {Wj^}  in  (62)  have  the 
common  exponential  PDF  used  in  example  (4)  for  all  the  earlier 
results  here.  Then,  by  reference  to  (31)  and  (39),  we  can  state 
that  the  common  PDF  of  RVs  fxj^l  defined  in  (62)  is 


-M-1 


p^(u)  =  [l  +  ^  for  u  >  0 


Similarly,  the  CDF  of  RVs  Ixj^)  follows,  by  reference  to  (38),  as 


23 


TR  8865 


C^(X)  =  1  -  [l  +  I]  ”  for  X  > 
The  CDF  of  RV  y  defined  in  (63)  is  then 


1  - 

1  Mj 

Cy(Y)  =  [C^(Y)]*^  = 


Its  inverse  function  follows  readily  as 

-1/M 


K 


for  Y  >  0 


-  1 


for  0  <  P  <  1 


(65) 


(66) 


Cy(P)  =  m[(i  -  p^/^) 

The  PDF  of  RV  y  can  be  obtained  by  differentiation  of  (66): 


(67) 


Py(u)  =  K  [C^(u)]^-1  p^(u)  = 


=  K 


1 


for  u  >  0  . 


(68) 


The  j-th  moment  of  RV  y  can  be  found  from  the  integral 

00 

Py(j)  =  J  du  u^  Py{u)  for  j  <  M  .  (69) 

0 

where  we  have  taken  advantage  of  the  fact  that  RV  y  can  never  be 
negative.  However,  a  useful  alternative  in  some  cases  is 
afforded  by  employing  integration  by  parts  on  (69),  with  the 
result,  for  j  2;  1,  that 
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Py(j)  =  j  J  du  ^  [1  -  Cy(u)]  ,  (70) 

0 

where  we  assume  that 

lim  u^[l  -  C  (u) ]  =  0  .  (71) 

u-^+®  ^ 

This  requirement  is  tantamount  to  presuming  that  the  j-th  moment 
Py(j)  exists,  that  is,  j  ^  M. 

When  we  use  CDF  (66)  for  RV  y,  then  j-th  moment  (70)  becomes 


00 

/Jy(j)  =  j  J  du  u^~^  < 

1  - 

1  ^  h) 

K' 

0 

In  order  to  evaluate  this  integral,  let,  for  the  moment, 

/  \  -M 

2  =  ^  S)  ■ 

Then,  the  bracketed  term  in  (72)  can  be  expanded  according  to 


k=0  k=l 


Employment  of  this  result  in  (72)  yields 
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=  jl 


M 


EEI 


Lzll 


k-1 


(Mk-j)j 


for  1  1  j  <  M  , 


(75) 


where  we  employed  [1;  3.194  3  and  8.384  1]  and  simplified  the  end 
result.  For  K  =  1,  (75)  reduces  to  (41),  as  it  must. 

The  first  two  moments  for  RV  y  follow  readily  from  (75)  as 

Py(l)  =  (*)  (76) 

k=l 

and 

K 

My(2)  =  2m2  0  (h^:2)(HK-1)  M  >  2  .  (77) 

k=l 


Both  of  these  results  are  easily  computed  by  means  of 
recurrences;  however,  a  bad  feature  of  these  two  sums  is  that 
they  are  alternating  series  and  lose  significance  for  large  K  due 
to  the  binomial  coefficient  which  gets  very  large.  A  backup 
procedure  is  to  revert  to  numerical  integration  of  (68)  -  (69)  or 
(70)  -  (72),  both  of  which  integrands  can  never  go  negative  and 
which  decay  as  u-'  for  large  u. 


26 


TR  8865 


However,  better  alternatives  to  the  first  and  second  moments 
have  been  found,  that  are  not  subject  to  cancellation  and  loss  of 
significance.  Namely,  it  is  shown  in  appendix  A  that 


Py(l)  =  M[F^  -  1]  , 


where 


for  M  >  2  . 


(78) 


(79) 


These  finite  products  are  very  useful  and  retain  significance 
even  for  large  K,  where  (76)  and  (77)  are  useless.  The  very 
compact  BASIC  program  listed  below  computes  both  quantities  in 
(78).  The  program  has  been  written  so  as  to  avoid  overflow, 
even  when  K  is  large. 


K= 

M=  !  M  >  2 

A=l/M 

B=2/M 

F1=F2=1 

FOR  Ks=l  TO  K 

Fl=Fl*Ks/(Ks-A) 

F2=F2*Ks/(Ks-B) 

NEXT  Ks 
Muy=M*(Fl-l) 

Sigy=M*SQR(F2-Fl*Fl )  (80) 


By  using  the  techniques  in  appendix  A  and  laboriously 
expanding  out  (75)  for  j  =3,  it  has  been  found  that 


/^y(3)  = 


M^(F. 


-  3F„  +  3F,  -  1) 


where  we  define  products 
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for  0  i  m  <  M  . 


Notice  that  =  1.  Based  on  this  result  and  (78),  we 
conjectured  that  the  j-th  moment  is  generally  given  by 


(n)  ^j-n  ior  j  <  H  . 

n=0 

In  fact,  this  has  been  confirmed  numerically  for  several  values 
of  M,  K,  and  j.  For  large  j,  the  alternating  character  of  this 
series  would  also  suffer  from  loss  of  significance;  however,  for 
the  low  order  moments  of  general  interest,  this  is  not  a  problem. 

The  third  and  fourth  cumulants  of  RV  y  were  also  computed  in 
terms  of  sequence  (F^j;  they  turned  out  to  be 

Xy(3)  =  M^[f3  -  3F2F3  +  2fJ]  , 

Xy(4)  =  M^[f4  -  -  3fI  +  12F2f2  -  6fJ]  . 

But  these  rules  for  obtaining  cumulants  {Xy(j)}  from  products 
{Fj^}  are  identical  to  the  general  rules  for  obtaining  cumulants 
from  moments,  within  the  factor  ;  see,  for  example, 

[4;  page  70,  (3.41)].  Thus,  we  have  a  very  efficient  and 
accurate  method  for  obtaining  the  low-order  cumulants  directly 
from  the  finite  products  The  case  for  Xy(l)  (78) 

differs  slightly  from  the  usual  rule,  in  the  need  to  subtract  1. 
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ANALYTICAL  CHECKS 

Numerous  checks  on  the  results  above  are  possible.  When 
M  =  «>  (no  normalization),  the  CDF  of  y  in  (66)  reduces  to  (57), 
where  it  must  be  remembered  that  K  in  this  section  on  combined 
normalization  and  or-ing  corresponds  to  M  in  the  section  on 
or-ing  alone.  On  the  other  hand,  if  K  =  1  (no  or-ing)  in  (66), 
the  result  in  (38)  correctly  emerges. 

With  regard  to  the  inverse  CDF  in  (67),  it  reduces  to  (58) 
for  M  =  00,  whereas  it  reduces  to  (40)  for  K  =  1.  The  PDF  of  RV 
y,  given  in  (68),  reduces  to  the  derivative  of  (57)  for  M  =  oo, 
whereas  it  reduces  to  (39)  for  K  =  1.  Finally,  the  first  two 
moments  in  (78)  -  (79)  reduce,  after  some  manipulations,  to  (56) 
for  M  =  oo;  on  the  other  hand,  the  j-th  moment  for  K  =  1  is  best 
handled  from  form  (75)  which  correctly  reduces  to  (41). 

EXTENSIONS 

The  case  where  the  normalizer  and  or-ing  device  are  followed 
by  an  averager  is  discussed  in  the  summary  below,  and  the  method 
of  determining  the  performance  of  that  system  is  outlined. 

Another  alternative  with  practical  merit  is  that  of  normalization 
followed  by  averaging  and  or-ing.  Since  the  CF  of  the  normalizer 
output  is  available  by  a  Fourier  transform  of  (39),  it  can  be 
raised  to  a  power  to  find  the  CF  at  the  averager  output.  Then 
another  Fourier  transform  can  yield  the  CDF.  At  this  point,  we 
could  utilize  (47)  to  find  the  CDF  of  the  system  output. 
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CHANGES  IN  SIZES  OF  NORMALIZER  AND  OR-ING  DEVICE 


We  now  address  the  change  in  size  of  the  normalizer  in  (62) 
from  size  M  to  N  and  the  change  in  size  of  the  or-ing  device  in 
(63)  from  K  to  L.  There  are  no  restrictions  on  the  relative 
sizes  of  the  parameters.  The  new  equations  are 


N  T'  '^kn  ^ 
n=l 


for  1  i  k  1  L  , 


(81) 


z  =  . .  ,yj^) 


The  CDF  of  RV  z  follows,  by  similarity  of  form  to  (66),  as 


C^(Z)  =  1  -  (l  +  I]  **  for  Z  >  0  . 


(83) 


The  inverse  function  is  easily  shown  to  be 


-1/N 


C^(P)  =  n[(i  -  pi/")  '  -  l] 


-  1  for  0  <  P  <  1 


The  first  two  moments  of  RV  z  in  (82)  are  given,  by 
comparison  with  (78)  and  (79),  as 


p^(l)  =  N[Gj  -  1]  ,  =  n(g2  -  gJ]  , 


(85) 


where 


1.1  IJ 

1 '  riT^l '  *^2 '  ft 


- 1  ■ 
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INPUT  SIGNAL-TO-NOISE  RATIO  REQUIREMENTS 

In  this  section,  we  will  determine  the  CDFs  of  the  outputs  of 
the  three  nonlinear  systems  considered  above,  namely  (31),  (46), 
and  (62)  -  (63),  for  the  case  where  a  signal  is  present  in  one  of 
the  input  RVs.  This  will  enable  us  to  use  the  results  given  in 
(27)  -  (29)  for  determination  of  required  input  SNRs  for  a 
specified  system  detection  probability  Pdj^  =  1-S^. 

NORMALIZER 


The  nonlinear  transformation  of  immediate  interest  is  the 
normalizcr  given  by  (31).  The  PDF  of  denominator  average  RV  a  is 
given  by  (37),  while  the  PDF  of  numerator  RV  Xq  with  signal 
present  will  be  modified  from  (4)  to  the  form 

^''pItTr)  for  u  >  C  ,  {87) 

where  R  is  the  input  power  SNR.  The  corresponding  CDF  is 

C^(X;R)  =  1  -  exp(^)  for  X  >  0  .  (88) 

By  reference  to  (33),  (38),  and  (88),  we  find  the  CDF  of  RV  y 
in  (31)  for  signal  present  to  be 


,M  ^M-1 


0 

=  ^  V  >  0  . 


(89) 
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We  now  employ  this  result  in  (27)  to  obtain 


S 


1 


Yt/M 


1+R. 


-M 


for  1  1  t  1  T  . 


(90) 


The  solution  for  the  required  input  SNR  R^  is  then  given  by 


Y  /M 

K  =  - TiTm -  -  1  for  1  i  t  1  T  .  (91) 

^  ( 1  -  Sj)  -  1 


We  must  repeat  here  the  caution  mentioned  in  (29),  namely 
that  specified  probability  at  threshold  Y^  must  be  less  than 
or  equal  to  probabilities  (P^l  in  order  that  nonnegative  SNRs 
{R^}  result  in  (91).  This  is  reasonable  since  it  corresponds 
physically  to  demanding  a  larger  detection  probability  when 
signal  is  present  than  when  absent. 
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OR-ING  DEVICE 


Here,  we  are  interested  in  the  or-ing  device  characterized  by 
(46)  when  signal  is  present  in  one  of  the  RVs  CDF  of 

RV  y  is  then  given  by 


Cy{Y;R)  =  C^(Y;R)  [C^(Y)  ]' 


=  1  -  exp 


[ 1  -  exp(-Y) ]  for  Y  >  0  , 


where  we  used  (88)  and  (4). 

If  we  now  employ  (92)  in  (27),  we  have  to  satisfy 


Sj  =  1  -  exp 


—  Y  '  1  M  “  1 

[l  -  exp(-Y^)]  for  1  <  t  <  T  .  (93) 

»  t  J  J 


The  solution  for  the  required  input  SNR  is  given  by 


^t  -  ln(l  -  Sj/Q^)  ^  for  1  1  t  1  T  , 


where 


2t  -  [l  - 


for  1  1  t  1  T 


Again,  (29)  must  be  satisfied. 


33 


TR  8865 


NORMALIZER  AND  OR-ING  DEVICE 

The  nonlinear  transformation  to  be  investigated  here  is  the 
combination  of  normalization  and  or-ing,  as  characterized  by  (62) 
and  (63),  when  signal  is  present  only  in  RV  Therefore,  only 

RV  in  (63)  contains  signal. 

The  CDFs  of  RVs  for  2  1  k  1  K,  are  given  by  (65).  On 

the  other  hand,  the  CDF  for  Xj^  is  available  by  reference  to  (89), 
in  the  form 


C^(X;R)  =  1  -  (l  +  for  X  >  0  . 


(96) 


Therefore,  the  CDF  for  RV  y  in  (63)  is  given  by 
Cy(Y;R)  =  C^(Y;R)  tCj^(Y)  1*^'^  = 


[l  -  fl  +  y/M]-"] 

1  -  [l  +  n] 

[  1+Rj 

[  mJ 

K-1 


for  Y  >  0  ,  (97) 


where  we  used  (96)  and  (65). 

When  we  equate  this  result  to  the  specified  probability 
according  to  (27),  we  obtain 


Y  /M' 

-M 

[  ^tl 

-M 

rv””  1. 

Si  = 

1  - 

^  1+R 

1  - 

1  ^  ir 

for  1  <  t  1  T  . 

tj 

k  * 

The  solution  for  the  required  input  SNRs  |R^1  is  given  by 
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R 


t 


(1  -  -  1 


1  for  1  i  t  1  T  , 


(99) 


where 


1  - 

-M 

K-1 


for  1  i  t  1  T 


(100) 


Restriction  (29)  must  be  satisfied  here  also. 

Finally,  if  several  detection  probabilities  Pd^^,  Pd2, 
are  of  interest,  we  must  satisfy  (30),  where 

max{l-P. 1  1  minlPd-l  .  (101) 

^  J 

As  checks  on  the  results  in  this  subsection  on  combined 
normalization  and  or-ing,  we  observe  for  M  =  »  (that  is,  no 
normalization),  (100)  reduces  to  (95),  where  K  here  corresponds 
to  M  there  for  or-ing  alone.  Also,  (99)  reduces  to  (94).  On  the 
other  hand,  for  finite  M,  but  with  K  =  1  (that  is,  no  or-ing), 
then  (100)  reduces  to  =  1,  in  which  case  (99)  reduces  to  (91). 
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SUMMARY 

We  have  determined  the  false  alarm  and  detection 
probabilities  for  three  different  nonlinear  transformations, 
namely  a  normalizer,  an  or-ing  device,  and  a  normalizer  followed 
by  or-ing.  In  particular,  results  are  given  for  the  following 
statistics  of  the  outputs  of  each  device:  the  PDF,  the  CDF,  the 
inverse  CDF,  and  either  the  moments  or  the  cumulants  (depending 
on  their  relative  tractability ) .  These  results  are  sufficient  to 
compute  the  receiver  operating  characteristics  (ROCs)  of  the 
processors,  if  desired.  However,  we  have  also  been  able  to  solve 
explicitly  for  the  input  SNR  required  to  realize  specified  false 
alaiir.  and  detection  probabilities;  this  largely  circumvents  the 
need  to  compute  ROCs. 

Two  programs,  with  numerical  examples  of  their  execution,  are 
listed  in  appendix  B.  The  first  corresponds  to  the  case  where 
the  number  T  of  thresholds  is  kept  fixed  as  the  size  of  the 
nonlinear  transformation  is  changed  from  M  to  N;  see  (20)  -  (23). 
On  the  other  hand,  the  second  program  allows  the  number  of 
thresholds  to  change  from  T  to  U  as  the  size  of  the  nonlinear 
transformation  is  changed  from  M  to  N;  see  (24)  -  (26). 

Both  programs  are  written  for  the  general  case  where  there  is 
both  normalization  (of  size  M)  and  or-ing  (of  size  K)  included  in 
the  data  processing;  see  (62)  -  (63).  By  making  M  infinite,  the 
program  will  handle  the  case  of  or-ing  alone;  on  the  other  hand, 
by  setting  K  =  1,  the  program  addresses  the  case  of  normalization 
alone.  Thus,  these  two  programs  cover  all  the  cases  addressed  in 
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this  investigation.  (Since  it  is  not  possible  to  actually  set 
M  infinite  in  a  program,  this  situation  is  handled  by  setting  M 
to  any  value  less  than  or  equal  to  2,  in  order  to  flag  this  case 
in  the  program,  and  then  branching  appropriately  at  various 
points  in  the  program.  The  substitute  equations  for  this  case  of 
infinite  M  come,  of  course,  from  the  earlier  analysis  for  or-ing 
alone.  For  finite  variance,  normalization  requires  M  >  2;  see, 
for  example,  (77)  or  (79).) 

We  have  not  included  the  effects  of  averaging  after  the 
normalization  and/or  or-ing  in  this  study.  Hence,  the  required 
input  SNRs  calculated  here  sometimes  turn  out  to  be  rather  large. 
The  exact  analysis  including  averaging  would  be  rather  involved, 
since  the  new  decision  variable  would  have  a  characteristic 
function  given  by  a  power  of  the  characteristic  function  f^il)  of 
current  output  variable  y;  that  is,  from  (68), 


(1) 


K  J  du  exp(i?;u) 
0 


1 


(102) 


This  is  probably  best  handled  through  the  use  of  fast  Fourier 
transforms.  The  integrand  decays  as  u~  ~  for  large  u,  which  is 
attractive  since  M,  the  normalization  size,  is  generally  fairly 
large  in  order  to  guarantee  decent  performance. 

In  the  meantime,  in  order  to  get  a  rather  rough  idea  of  the 
improvement  attainable  by  employment  of  averaging,  it  is 
suggested  that  the  rule  of  thumb  (3j  (C-10)]  for  the  input 
signal-to-noise  ratio  improvement  in  dB,  -5  log  A,  be  used,  where 
A  is  the  number  of  independent  quantities  averaged. 
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APPENDIX  A.  SIMPLIFICATION  OF  SUMS  (76)  AND  (77) 

Here,  we  will  convert  the  alternating  series  in  (76)  and  (77) 
into  finite  products  that  retain  significance,  even  for  large 
values  of  K.  We  begin  with  the  first  line  of  (52)  and  expand  the 
power  term  in  a  binomial  series,  obtaining 

00  M-1 

=  M  J  du  exp(ij;u-u)  (-1)^  exp(-ku)  = 

0  feo 


M-1 


(-1) 


k  +  1  -  i^ 


(A-1) 


Now,  equate  (A-1)  to  its  alternative  expression  in  the  last 
line  of  (52),  and  then  replace  M  everywhere  by  K+1,  getting 


(K+1) 


EE) 


k 


+ 


(A-2) 


Now  let  z  =  1  -  iJ;  in  (A-2)  and  simplify;  there  follows 


HIk) 


(-1) 
k  +  z 


K! 


(z) 


K+1 


{A-3) 


Thus,  the  alternating  series  has  been  converted  into  a  finite 
product  which  is  useful  for  all  values  of  K  without  loss  of 
significance. 
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The  use  of  (A-3)  on  (76)  yields  the  result  quoted  in  (78)  and 
(79)  for  the  first  moment  /jy(l).  On  the  other  hand,  in  order  to 
convert  (77),  it  is  necessary  to  take  the  preliminary  step  of 
breaking  the  rational  function  into  two  parts  according  to 


1  ^  [  JL_  1  1 

(k-b)(k-a)  b-a  Ik-b  "  k-aj  " 


(A-4) 


and  then  to  use  (A-3)  once  with  z  =  -1/M  and  once  with  z  =  -2/M. 
After  manipulation,  simplification,  and  cancellation  of  common 
terms  with  the  square  of  Py(l),  the  end  result  for  the  standard 
deviation  of  RV  y  is  found  to  be  just  the  second  term  in  (78). 

The  results  in  (78)  -  (79)  have  been  numerically  checked  with  the 
original  defining  integral  (72),  for  j  =  1  and  j  =  2,  for  several 
values  of  M  and  K;  they  agree  within  the  accuracy  of  the  computer 
employed. 

A  more  general  result  than  (A-3)  is  available  by  means  of  a 
different  approach.  First,  for  general  values  of  a,  note  that 


al  ^  (-1)*^  ^~^)k  1  ^  1  ^^^k 

kj  kl  '  k+z  z  (z+l)j^ 


(A-5) 


Then,  the  following  alternating  sum  can  be  manipulated  into  a 
more  useful  form,  namely 


(-1)*"  =  1 
k+z  z 


00 


(-a)k  (2)^ 

(z+1)^  k! 


^  F(-a,z;z+l; 1 ) 


r(l+a)  r(z) 
r(l+a+z)  ' 

(A-6) 


where  we  used  [2;  15.1.1  and  15.1.20].  If  we  set  a  =  K  in 
(A-6),  it  reduces  to  (A-3).  (A-6)  is  also  equal  to  B(l+a,z). 
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APPENDIX  B.  PROGRAMS 

There  a-O  two  frojraiP'-  in  thin  appendix,  both  in  BASIC 

for  the  Hewlett-Packard  9000  computer.  None  of  the  variables  are 
declared  INTEGER;  thus,  for  example,  input  parameters  M,  K,  N,  L 
are  all  treated  as  REAL  variables  throughout. 

The  first  program,  listed  on  pages  42  -  44,  requires  that  the 
number  of  thresholds  T  be  maintained  the  same  when  the  sizes,  M 
and  K,  of  the  normalizer  and  or-ing  device,  respectively,  are 
changed  to  N  and  L.  On  the  other  hand,  the  second  program, 
listed  on  pages  46  -  49,  allows  the  number  of  thresholds  to 
change  from  T  to  U,  which  can  be  either  larger  or  smaller. 

The  listings  are  heavily  keyed  to  the  explicit  equations  in 
the  main  text;  this  should  enable  the  user  to  identify  and  modify 
particular  manipulations  if  desired.  It  should  be  noted  that,  in 
the  programs,  the  or-ing  size  begins  at  value  K  and  gets  changed 
to  L.  If  these  results  are  to  be  compared  with  the  or-ing  only 
results  in  the  text,  where  the  parameter  M  was  used,  it  is 
necessary  to  make  the  switch  from  K  in  the  program  to  M  in  the 
text.  Example  outputs  from  both  programs  are  presented  after  the 
listing  for  each  case. 
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10 

!  NORMFlLIZER  ORING,  EQUISPRCED  IN  POWER,  SRME  NO.  OF  THRESHOLD 

20 

Pl=. 85 

j 

SPECIFIED  PROBABILITY,  (9) 

30 

Dby=10 

1 

DECIBEL  RATIO  OF  THRESHOLDS,  <10> 

40 

T  =  7 

\ 

HUMBER  OF  THRESHOLDS,  (10) 

50 

M  -  8 

\ 

INITIAL  NORMAL IZER  SIZE  >  2,  <62) 

60 

!  FOR  NO  INITIAL  NORMALIZATION,  THAT  IS,  M  INFINITE,  SET  M  <=  2 

70 

K  =  1  2 

! 

INITIAL  OR-ING  SIZE  >  0,  <63) 

80 

N=16 

1 

NEW  NORMRLIZER  SIZE  >  2,  <81 ) 

90 

!  FOR  HO  NEW  NORMALIZATION, 

THAT  IS,  N  INFINITE,  SET  N  <=  2 

1  00 

L  =  24 

1 

NEW  OR-ING  SIZE  >  0,  <82) 

1  10 

Pdl=. 5 

I 

SPECIFIED 

1  20 

Pd2=. 9 

1 

DETECTION 

130 

Pd3=. 99 

1 

PROBABILITIES,  <30) 

140 

PRINT  "PI  =";P1;" 

Dby 

=";Dbyj"  T  =";T 

150 

PRINT  "M  =";N;" 

K 

_1 

II 

-J 

II 

II 

160 

REDIM  Y<:i;T),Yb<l:T>, 

P<  1 : 

T),2<l:T),Zb<l:T) 

170 

DIM  Y<99)  ,  Yb<99) , P<99> , Z<99> , 2b<99> 

180 

Rk  =  l.'K 

1 

<67) 

190 

R1 =1/L 

! 

<84) 

200 

IF  M>2  THEN  300 

210 

F1=F2  =  0 

! 

M  INFINITE 

220 

FOR  Ms=l  TO  K 

1 

<56),  m 

230 

Rl  =  l.  Ms 

240 

F1=F1+R1 

250 

F2=F2+R1*R1 

260 

NEXT  Ns 

270 

Muy  =  F  1 

! 

<56) 

280 

S 1 gy  =  SQR  <  F2 ) 

! 

<56> 

290 

GOTO  390 

300 

Rl=l/M 

! 

M  >  2 

310 

R2  =  2-'M 

320 

F1=F2=1 

330 

FOR  Ks=l  TO  K 

1 

<79>,  k 

340 

Fl=Fl*Ks/-<Ks-Rl) 

! 

<79) 

350 

F2=F2*Ks'<Ks-R2) 

! 

<79) 

360 

NEXT  Ks 

370 

Nuy=N*<Fl-l) 

! 

<78) 

380 

Sigy=M*SQR<F2-Fl*Fl) 

! 

<78) 

390 

IF  W>2  THEN  490 

400 

G 1 =G2=0 

1 

N  INFINITE 

4  1  0 

FOR  Ns=l  TO  L 

! 

<61),  n 

420 

Rl=l/Ns 

430 

G1=GI+R1 

440 

G2=G2+R1*R1 

450 

NEXT  Ns 

460 

Mu2=G1 

! 

<61  ) 

470 

S i gz  =  SQR  <  G2  > 

! 

<61  ) 

480 

GOTO  580 

490 

Rl=l/N 

1 

N  >  2 

500 

R2=2/N 

510 

G1=G2=1 

520 

FOR  Ks=l  TO  L 

1 

<86),  k 

530 

G1=G1*Ks.^<:Ks-R1  ) 

1 

<86) 

540 

G  2  =  G  2  *  K  s  /  C  K  s  -  R  2  ) 

1 

<86) 

550 

NEXT  Ks 

560 

Muz  =  N*  <  G 1 -  1 ) 

t 

<85) 

570 

SI  g2  =  N*SQR(’G2-Gl*Gl  ) 

I 

<85) 

530 

PRINT  "Muy  =";ML4y;" 

Sigy  =";Sigy 

590 

PRINT  "Muz  =";Nuz5" 

Sigz  ="5Sigz 

600 

PRINT 
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610 

j 

EGlUISPflCEI)  IN  DECIBELS: 

REMOVE  880-970  AND  INSERT  620-710 

620 

1 

R=l-Pl''Rk  ! 

<9)  8.  <67) 

630 

! 

IF  ri>2  THEN  660 

640 

! 

Y1=-L0G<R)  ! 

(58) 

650 

1 

GOTO  670 

660 

1 

Y  ?  'M?  -  1  '  ! 

CS7) 

670 

! 

D1  =  10*LGT<Y1  )  1 

(14) 

680 

! 

Del  d  =  Dbi^/<T-l )  1 

(15)  &  (16) 

690 

1 

FOR  Ts=l  TO  T  1 

(16),  t 

700 

j 

Dt=Dl+Deld*<Ts-l )  ! 

(16) 

710 

! 

Y<Ts)=Y=10''(:Dt/10)  ! 

(17) 

720 

! 

PRESET  CONSTANTS:  REMOVE 

880-980  AND  INSERT  730-760 

730 

! 

DATA  1,3,5,7,9,10,11  ! 

USER  MUST  INPUT  T  NUMBERS 

740 

! 

READ  Ybf*)  I 

NORMALIZED  Y  THRESHOLDS 

750 

1 

FOR  Ts=l  TO  T  ! 

(18),  t 

760 

I 

Y<Ts)=Y  =  Muy  +  Si  .3y*Yb<Ts)  ! 

(18) 

770 

) 

PRESET  PROBABILITIES:  REMOVE  880-1030  AND  INSERT  780-870 

780 

1 

DATA  . 9, . 99, . 999, . 9999, . 

99999, . 999999, . 9999999 

790 

1 

READ  P<*)  1 

(19) 

800 

j 

FOR  Ts=l  TO  T  ! 

(19),  t 

810 

! 

P=P<Ts) 

820 

! 

R=l-P^Rk  ! 

(58) 

830 

1 

IF  N>2  THEN  860 

840 

1 

Y(:Ts)=Y  =  -LOG<R)  ! 

(58) 

850 

i 

GOTO  870 

860 

! 

Y<Ts)=Y=N*<R^<-l/M)-l )  ! 

(9)  8.  (67) 

870 

! 

Yb<T£)=<Y-Nuy)/Si gy  ! 

(12) 

880 

R=l-Pl'^Rk  ( 

(9)  8<  (67) 

890 

IF  M>2  THEN  920 

900 

Yl=-LOG<R)  ! 

(58) 

910 

GOTO  930 

920 

Y1=M*<R^<-1/M>-1 )  ! 

(67) 

930 

Ay=10''<:Dby/10)  1 

(10) 

940 

Ytc=Yl*Ay  ! 

(10) 

950 

Del y=<Ytc-Yl )/<T-l )  ! 

(11) 

960 

FOR  T£=1  to  T  I 

(11),  t 

970 

Y<T£)=Y=Y1+Del y*<Ts-l )  ! 

(ID 

980 

Yb<Ts)=<Y-Muy)/Si gy  ! 

(12) 

990 

IF  M>2  THEN  1020 

1000 

R=EXP<-Y)  ! 

(57) 

1010 

GOTO  1030 

1020 

R=<  1+Y/'M)''<-M)  I 

(13)  8.  (66) 

1030 

P<Ts)=P=<  1-R)''K  1 

(66) 

1040 

Q=1-P^R1  1 

(84) 

1050 

IF  N>2  THEN  1080 

1060 

2<T£)=2=-L0G<Q)  ! 

(21)  8.  (60) 

1070 

GOTO  1090 

1080 

2<T£)=2  =  N*<Q''<-1/N)-1  )  ! 

(21)  8.  (84) 

1090 

Zb < Ts )=< 2-Muz ) /S i gz  ! 

(22) 

1100 

NEXT  Ts 

1110 

Dbz=10»LGT<Z<:T)/Z<  1  )  )  1 

(23) 

1  1  20 

PRINT  "  Y  THRESHOLDS 

NORMALIZED  PROBABILITIES 

1  1  30 

FOR  T£=1  to  T 

1  140 

PRINT  Y(T£) , Yb<T£) , P<T£) 

1  1  50 

NEXT  Ts 

1160 

PRINT 

1  1  70 

PRINT  "  2  THRESHOLDS 

NORMAL  I  ZED  PROBABILITIES 

1  130 

FOR  T£=1  to  T 

1  190 

PRINT  2<:T£),2b(T£),P<Ts> 

1  200 

NEXT  Ts 

1210 

PRINT 

1220 

PRINT  "Dbz  ="jDbz 

43 


TR  8865 


1230 

PRIHT 

1240 

PRINT  "  Pdl  ='‘;Pdl;‘' 

Pd2  ="  ;  Pd2;  " 

Pd3  =";Pd3 

1250 

PRINT 

1260 

PRINT  "SIGNFlL-TO-NOISE 

RATIOS  <DB>  FOR  INITIAL  TRANSFORMATION:" 

1270 

IF  N>2  THEN  1400 

1280 

FOR  Ts=l  TO  T 

! 

<94),  t 

1290 

Y=Y<Ts) 

1300 

Q=(  l-EXP(-Y)  )^<:k-i  ) 

I 

<95) 

1310 

D1=-L0G< 1-C 1-Pdl 

I 

<94> 

1320 

Rt  1=Y''D1-1 

1 

<94) 

1330 

D2  =  -L0G(1-C  1-Pd2)/'Q) 

1340 

Rl2  =  Y-'D2-l 

1350 

D3=-L0G< l-( 1-Pd3)/Q) 

1360 

Rt3  =  Y.-D3-l 

1370 

PRINT  10*LGT<Rt.  1  >  ,  10»LGT<Rt2),  10*LGT<Rt.3) 

1380 

NEXT  Ts 

1390 

GOTO  1530 

1400 

Rtn  =  -1/M 

1 

(99) 

1410 

FOR  T£=1  to  T 

! 

(99),  t 

1420 

Ym  =  Y(Ts).^M 

1 

(99)  &  (100) 

1430 

Q=<;  1+Ytn  -M> 

! 

(100) 

1440 

Q=< l-Q^-CK-l ) 

! 

(100) 

1450 

Dl  =  a-<  1-Pdl  )/Q)^Rm-l 

1 

(99) 

1460 

Rt  l=Yfii/Dl-l 

! 

(99) 

1470 

D2=< l-C l-Pd2)/Q)^Rm-l 

1480 

Rl2=Yrrt/D2-l 

1490 

D3=<  l-(  l-Pd3)-  Q)''Rm-l 

1500 

Rt  3  =  Ym/  D3- 1 

1510 

PRINT  10*LGT<Rt 1 > , 10*LGT<Rt2), 10*LGT<Rt3) 

1520 

NEXT  Ts 

1530 

PRINT 

1540 

PRINT  "SIGNftL-TO-NOISE 

RATIOS  CDB)  FOR  NEW 

TRANSFORMATION: " 

1550 

IF  N>2  THEN  1680 

1560 

FOR  T£=1  to  T 

1 

SIMILAR 

1570 

2=Z<Ts) 

1 

TO 

1580 

Q=c:  l-EXPC-Z)  )^<L-1  ) 

1 

(94) 

1590 

D1=-L0G< l-< 1-Pdl )/Q) 

! 

& 

1  600 

Rt 1=Z/D1-1 

1 

(95) 

1610 

D2  =  -L0G<:  l-<  1-Pd2)/Q) 

1620 

Rt2=Z/D2-l 

1630 

D3  =  -L0GU-<  1-Pd3)/Q) 

1640 

Rt3  =  Z.'  D3-1 

1650 

PRINT  10*LGT<Rt 1 ) , 10*LGT<Rt2), 10*LGT(Rt3) 

1660 

NEXT  Ts 

1670 

GOTO  1810 

1680 

Rn=-1/N 

t 

SIMILAR 

1690 

FOR  Ts=l  TO  T 

! 

TO 

1  700 

Zn  =  Z';Ts)/N 

! 

(99) 

1710 

Q=(  l+2n>''<-N) 

1 

i, 

1720 

Q=(  l-CO-CL-l  ) 

1 

(100) 

1  730 

Dl  =  <  l-<  1-Pdl  )/^Q)''Rn-l 

1740 

Rt l=Zn/Dl-l 

1750 

I)2=(  l-<  l-Pd2)/Q)^Rn-l 

1  760 

Rt2  =  Zfv'D2-l 

1770 

D3=(:  l-a-Pd3:)/Q>-Rri-l 

1  780 

R  t  3  =  Z  n  z  D  3  -  1 

1  790 

PRINT  10*LGT'.'Rt  1  )  ,  1  0*LGT  <  Rt  2  )  ,  1  0«LGT  <  Rt  3  > 

1800 

NEXT  Ts 

1810 

PRINT 

1820 

END 
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=  .wt>  Dby  = 

10  T  =  7 

M  8  K 

=  12 

N  =  16 

L  =  24 

Nuy  =  3.94631506063 

Sigy  = 

2.09538979024 

Huz  =  4.32438211 

761 

S  i  gz  = 

1.69328544323 

Y  THRESHOLDS 

NORMALIZED 

PROBABILITIES 

5. 7085601983 

. 841010653879 

.  85 

14.2714004957 

4.92752493268 

. 996679075641 

22.8342407932 

9.01403921 147 

. 999753631023 

3 1 . 39708 1 0906 

13. 1005534903 

.99996531 1792 

39.9599213881 

17. 1870677691 

. 99999280763 

48.5227616855 

21 . 2735820479 

. 999998067526 

57. 085601983 

25. 360096326? 

. 999999374796 

2  THRESHOLDS 

NORMALIZED 

PROBABILITIES 

5. 66721821227 

.911 149446675 

.  85 

1 1 . 8779764405 

4.46091021044 

. 996679075641 

16.8024176865 

7. 3'=;9 12705339 

.999753631023 

.^.1 . 0784 9 26281 

9 . 894404  1 

. 999965311792 

24. 909916314 

12.  1571553566 

.  99999280763 

28.4120734509 

14.2254168838 

. 999998067526 

31 . 6575625522 

16.  1420985125 

. 999999374796 

Dbz  =  7. 32045232969 

Pdl  =  .5 

Pd2  =  .9 

Pd3  =  .99 

SIGNRL-TO-HOISE 

RATIOS  <DB)  FOR  INI 

TIAL  TRANSFORMATION 

7. 18029480646 

16.5239636654 

26.8808304974 

12.6997498763 

21 . 2426978346 

31 . 503743108 

1 4 . 8466399857 

23. 3091924052 

33.5584316626 

16.269595645 

24.6936419313 

34.9428284925 

17. 3390310658 

25.7492651765 

35.9905902914 

18. 1963923488 

26. 5945353622 

36.8340131363 

1 8 . 9120679996 

27 . 301 778399 1 

37. 5399641554 

SIGHHL-TO-HOISE 

RATIOS  <DE)  FOR  NEW 

TRANSFORMATION: 

7. 40601673703 

16. 6306689407 

26.9721139815 

1 1 . 95598 11212 

20.4527339701 

30. 7066466017 

13. 5624054665 

21 . 9841310487 

32.2269170012 

14. 587035072 

22.975309045 

33.21299621  17 

15.3357177535 

23.704135924 

33.9387696157 

15. 922698397 

24.2777465869 

34.5103143719 

16, 4  0  y  2  1  7  3  y  6 

24 . 749160188 

34.9802228898 
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10 
20 
30 
•4  0 
50 
60 
70 
80 
80 
100 
1  10 
1  20 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 


HORMRL I 
Pl=. 85 
Dbv‘=  1  0 
T  =  7 
U=15 
ri=8 

FOR  no 
k:=i2 

N=16 
FOR  NO 
L  =  24 
Pdl=. 5 
Pd2=. 9 
Pd3= . 99 
PRINT  "PI 


hK 


SPECIFIED  _ _ _  ,  . 

DECIBEL  RATIO  OF  THRESHOLDS,  ( 
INITIAL  NUMBER  OF  THRESHOLDS, 

NEW  NUMBER  OF  THRESHOLDS,  (24> 
INITIAL  NORMAL  I ZER  SIZE  >  2,  <62 ) 


j 

NEW  NORMALIZATION, 

I 


rh  iHi  "t-'i  =";P1;  ■■  lid^>  =";i.iDy;"  i  =" 

PRINT  "M  =";M;"  K  =";K;“  N  =";N;" 

RED  I M  V  a  :  T  > ,  Yb  <  1 !  T  i ,  P  <  I ;  T  )  ,  Zp  <  1 :  IJ  >  ,  Zbp  <  1 :  U  > ,  Pp  <  1 :  U ) 
D  I  M  Y  <  99  > ,  Yb  (  99  )  ,  P  <  99  ) ,  Zp<  99  )  ,  Zbp  <  99  >  ,  Pp  <  99  :> 

Rk=l^K  !  <67) 

Rl=l/L  !  (84) 


2, 

SET 
(63) 
(31  ) 
N  <  = 


NEW  NUMBER  OF  THRESHOLDS,  (24) 
!  INITIAL  NORMAL  I ZER  SIZE  >  2,  ( 
INITIAL  NORMALIZATION,  THAT  IS,  M  INFINITE,  SET  M 
!  INITIAL  OR-ING  SIZE  >  0,  (63) 

!  NEW  NORMAL  I ZER  SIZE  >  2, 

THAT  IS,  N  INFINITE,  SET 
NEW  OR-ING  SIZE  >  0,  (82 

SPECIFIED 
DETECTION 

PROBABILITIES,  (30) 

=  "  ;  Dbs^;  "  T  =" ;  T;  " 

;  N;  '• 


(10) 


U  = 
L 


D  I  M 

Rk  =  l^t( 

R1 =1/L 
IF  M>2  THEN 
F1=F2=0 
FOR  Ms=l  TO 


310 


K 


(84) 

M  INFINITE 
(56) ,  m 


240 

Rl  =  l/Hs 

250 

F1=F1+R1 

260 

F2=F2+R1*R1 

270 

NEXT  Ms 

280 

Mo(y  =  Fl 

!  (56) 

290 

Si gv-SQR(F2) 

!  (56) 

300 

GOTO  400 

310 

R1  =  1.''N 

!  M  >  2 

320 

R2  =  2/H 

330 

F1=F2=1 

340 

FOR  Ks=l  TO  K 

!  (79),  k 

350 

Fl=Fl*Ks/(Ks-Rl) 

!  (79) 

360 

F2=F2*Ks'(Ks-R2) 

!  (79) 

370 

NEXT  Ks 

330 

Muv=M*(Fl-l ) 

!  (78) 

390 

Si  gv  =  M*SG!R(F2-Fl*Fl  ) 

!  (78) 

400 

IF  N>2  THEN  500 

4  1  0 

G1  =G2  =  0 

!  N  INFINITE 

420 

FOR  Ns=l  TO  L 

!  (61),  n 

4  30 

Rl  =  l-'Ns 

440 

G 1 =G 1 +R 1 

450 

G2=G2+R 1 *R 1 

460 

NEXT  Ns 

470 

Muz=G 1 

!  (61  ) 

480 

Si gz=SQR(G2) 

!  (61 ) 

490 

GOTO  590 

500 

R1  =  1.-'N 

i  N  >  2 

510 

R2  =  2-'N 

520 

G1=G2=1 

530 

FOR  Ks=l  TO  L 

!  (86),  k 

540 

Gl=Gl*Ks^(Ks-Rl  ) 

1  (86) 

550 

G  2  =  G  2  •*  K  s  (  K  s  -  R  2  ) 

!  (86) 

560 

NEXT  Ki 

5  70 

M  Li  z  =  N  *  (  G  1  -  1  ) 

!  (85) 

580 

Si  g  z  =  N  *  S  0  R  (  G  2  -  G  1  *  G  1  ) 

1  (85) 

590 

PRINT  "Muy  =";Muy;" 

Sigy  ="5 

600 

PRINT  "Mvtz  =";Muz;" 

Sigz  ="; 

6  1  0 

PRINT 
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bZt)  !  EQUli.PftCED  IH  DECIBELS;  REMOVE  890-980  RMD  INSERT  630-728 
630  !  R=l-Pl''Rk  !  <9)  8,  <67> 

640  !  IF  M>2  THEN  670 

650  !  Y1=-L0G<:R)  !  (58) 

660  !  GOTO  680 

670  !  Yl=M*<R-(-l/M)-l )  !  (67) 

680  !  D1  =  10*LGT(  Y1 )  !  (14) 

690  !  De 1 d=Dby>( T- 1 )  !  (15)  &  (16) 

700  !  FOR  Ts=l  TO  T  t  (16),  t 

710  !  Dt =D1+De1 d*(T£-l )  !  (16) 

720  !  Y(T£)=Y=10' (Dt/10)  !  (17) 

730  !  PRESET  CONSTANTS;  REMOVE  890-990  AND  INSERT  740-770 

740  !  DATA  1,3,5,7,9,10,11  !  USER  MUST  INPUT  T  NUMBERS 

750  !  READ  Yb(*)  !  NORMALIZED  Y  THRESHOLDS 

760  !  FOR  Ts=l  TO  T  I  (18),  t- 

770  !  Y(T£)=V=Nuy+Si gy*Yb(T£) !  (18) 

780  !  PRESET  PROBABILITIES:  REMOVE  890-1040  AND  INSERT  790-8S0 

790  !  DATA  .  9 ,. 99 ,. 999 ,. 9999 ,. 99999 ,. 999999 ,. 9999999 

800  !  READ  P(*)  !  (19) 

810  !  FOR  T£=1  to  T  !  US),  t 

820  !  P-P(T3) 

830  !  Rsl-P'^Rk  !  (58) 

840  !  IF  N>2  THEN  870 

850  !  Y(Ts)=Y=-LOG(R)  !  (58) 

860  !  GOTO  880 

870  '  Y(T£)=Y  =  M»(R-^(-1.'N)-1  )  !  (9)  8.  (67) 

880  !  Yb(T£)=(Y-Nuy)/Si gy  I  (12) 

890  R=l-Pl-Rk  !  (9)  8<  (67) 

900  IF  H>2  THEN  930 

910  Y1=-L0G(R)  !  (58) 

920  GOTO  940 

930  Y1=M»(R'^'(-1/M)-1  )  I  (67) 

940  Ay=10''(Dby/10)  I  (10) 

950  Ytc=Yl*Ay  !  (10) 

960  Del y=(Yic-Yl )/(T-l )  !  (11) 

970  FOR  T£=1  to  T  !  (11),  t 

980  Y(T3)=Y  =  Y1+De1 y*(T£-l )  !  (11) 

990  Yb(Ts)=(Y-Muy)/Si gy  !  (12) 

1000  IF  M>2  THEN  1030 

1010  R=EXP(-Y)  !  (57) 

1020  GOTO  1040 

1030  R=(  1+Y^'M>''(-H)  !  (13)  8.  (66) 

1040  P(T£)=( 1-R)'K  !  (66) 

1050  NEXT  Ti 
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1  060 

Q1=1-P(1)'R1  ! 

(84) 

1070 

Of  c  =  l-P<T)  ■'Rl  ! 

(84) 

1  080 

IF  N : 2  THEN  1120 

1090 

Z  1  =  -  L  0  G  v  Q  1  )  ! 

(60) 

1  100 

Z  t  c  =  -  L  0  G  (  0 1  c  )  ! 

(60) 

1110 

GOTO  1150 

1  1  20 

Rn=-1/N  ! 

(84) 

1  1  30 

2l=N*cQr  Rn-1  )  > 

(21  ) 

& 

(84) 

1  140 

Ztc=N*u9tc  Rn-1)  ! 

(21  ) 

& 

(84) 

1  150 

Del  z  =  (Zt  c-Zl  )/-  (U-l  )  ! 

(24) 

1  160 

FOR  Us=l  TO  U  ! 

(24), 

u 

1  170 

Zp':  Us  )=Z  =  Zl+Del  z*(Lls-l  )  ! 

(24) 

1  1  SO 

Z  b  p  LI  s  )  =  <  Z  -  M  1.4  z  )  -'  3  i  g  z  ! 

(26) 

1  190 

IF  N>2  THEN  1220 

1200 

Pp<Us)  =  <:  l-EXP(-Z)  )  -  L  ! 

(25) 

(59) 

1210  GOTO  1230 

1220  PpOJi  >  =  (  l-U+2/N)''<-fO)-'L  !  (25)  t.  (83) 

1230  NEXT  Us 

1240  ribz=10*LGT(2p(U)/2p(  1  ))  !  (23) 

1250  PRINT  "  Y  THRESHOLDS  NORMALIZED  PROEhE  I L I T I ES '' 

1260  FOR  Ts=l  Til  T 

1270  PRINT  Y( Ts ) , Yb(Ts) , P(Ts) 

1280  NEXT  Ts 

1290  PRINT 

1300  PRINT  "  Z  THRESHOLDS  NORMALIZED  PROEAE I L I T I ES " 

1310  FOP  Us=l  TO  U 

1320  PRINT  ZpdJs) , 2bp(Us) , Pp(Us) 

1  330  NEXT  Us 

1340  PRINT 

1350  PRINT  "Dbz  =";Dbz 

1 3c.O  PRINT 

1 3  fa  1  PAUSE 

1  370  PRINT  "  Pdl  =";Pdl!‘'  Pd2  =";Pd2;"  Pd  3  ='';Pd3 

1380  PRINT 

1390  PRINT  "SIGNAL-TO-NOISE  RATIOS  (DB)  FOR  INITIAL  TRANSFORMATION:" 
1400  IF  rr;2  THEN  1530 

1410  FOR  Ts=:l  TO  T  !  (94),  t 

14  20  Y  =  ■(  (  T  s  ) 

1430  Q=':  l-EXPil-Y)  )^(K-1  )  !  (95) 

1440  D1=-L0G( l-< 1-Pdl )/Q)  !  (94) 

1450  Rtl=YVDl-l  !  (94) 

1460  D2=-L0Gc l-( 1-Pd2)/Q) 

14  70  R  t.  2  =  Y  /  D  2  -  1 

14  80  D  3  =  - L  0  G ( 1  - < 1  - P  d  3 ) / Q ) 

14  90  R  t  3  =  Y / D  3  -  1 

1500  PRINT  10*LGT' RM ) , 10*LGT(Rt2) , 10»LGT(Rt 3) 

1510  NEXT  Ts 

1520  GOTO  IbfaO 
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1  530 
'  540 
1550 
1560 
1570 
1580 
1590 
1  600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1  680 
1690 
1700 
1710 
172C 
1730 
1740 
1750 
1760 
1770 
1  780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
1  920 
1930 
1  940 
1950 


Rrii  =  -  1  /M 

1 

<99> 

FOR  Ts=l  TO  T 

j 

<99),  t 

Y(i.  =  Y':.Ts>/H 

1 

<99)  e,  (100> 

Q=<  1+Yni> '  C-M) 

1 

(100) 

Q=< 1-Q>^<K-1  ) 

! 

ClOO) 

Dl  =  <  l-<  1-Pdl  )/Q)''Rt(.-l 

1 

<99> 

Rt  l=Ym/'Dl-l 

D2=>:  l-<  1-Pd2)-'Q)-Rrn-1 
Rt2=Ym/D2-l 

D3=i'  l-<  1-Pd3)/Q)-Rri)-1 

! 

<99) 

Rt  3  =  Ym/D3-l 

PRINT  10*LGT<Rt  1  >  ,  10*LGT<Rl2)  ,  10»LGT(:Rl3> 

NEXT  Ts 
PRINT 

PRINT  "SIGNFIL-TO-NOISE  RRTIOS  (DB>  FOR  NEW  TRFiNSFuRMFiT  I  ON ;  " 


IF  H>2  THEN  1810 

FOR  Us=l  TO  U 

!  SIMILAR 

Z=2p<Us) 

!  TO 

Q=< 1-EXP<-Z) J-CL-l  > 

!  <94) 

D1=-L0G<  l-a-Pdl)/Q) 

!  t 

R1=2-'D1-1 

!  (95) 

D2=-L0G< l-< 1-Pd2)/Q) 
R2=Z/D2-1 

D3  =  -L0G<  l-<  1-Pd3)/G!) 
R3=Z/D3-1 

PRINT  10*LGT<R1), 10*LGT<R2), 10*LGT(R3) 


NEXT  Us 
GOTO  1940 

Rn  =  -1/'W  !  SiniLflR 

FOR  Us=l  TO  U  !  TO 

2ri  =  Zp<Us)/N  !  <99) 

Q=<  1+Zn)''<-N)  !  & 

Q=<  1-Q>''(L-1)  !  aOO) 


Dl  =  <  l-<  1-Pdl  )/Q)''Rri-l 
Rl=Zri/Dl-l 

D2=<  l-<  l-Pd2)/Q)^Rri-l 
R2=Zn/D2-l 

D3=(  l-<  l-Pd3)/Q)''Rn-l 
R3  =  Zri/D3-  1 

PRINT  10*LGT<R1 ) , 10*LGT<R2) , 10*LGT<R3) 

NEXT  Us 

PRINT 

END 
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P  1  =  .  3  5  Ii  b  y  = 

10  T  =  ? 

U  =  15 

M  =  8  K  =  12 

N  =  16 

L  =  24 

Muy  =  3.84831506068 

Sigy  = 

2.09 

538879024 

Muz  =  4.32438211761 

S  i  gz  = 

1 . 69328544323 

Y  THRESHOLDS 

NORMALIZED 

PROBABILITIES 

5. 7085601883 

.  34  1010653878 

.  85 

14. 2714004857 

4.82752483268 

. 996678075641 

22.8342407832 

8.01403821147 

. 999753631023 

3 1 . 3870810806 

1 3. 1005534803 

.999965311782 

38 . 85882 1 3881 

17. 1870677681 

. 99899280763 

48.5227616855 

21 . 2735820478 

.998883067526 

57. 085601883 

25. 3600863267 

. 999998374796 

2  THRESHOLDS 

NORMALIZED 

PROBABILITIES 

5.86721821227 

.81 1148446675 

.  85 

7 . 708385665 1 2 

1 . 88807437885 

. 956528647638 

8.55155311787 

3.08688831323 

.98667504332 

11.5837205708 

4. 1748242465 

. 99560662072 

1  3 . 2  3  5  3  8  8  0  li  3  7 

5. 26284817878 

. 998447087585 

15.0780554765 

6.35077411306 

.899415548265 

16. 8202228284 

7. 43868804633 

.999767363279 

18. 7623803822 

8. 52662387861 

. 999902645376 

20.6045578351 

8.61454881288 

.999957385298 

22. 4467252878 

10.7024738462 

. 999980574274 

24.2888827408 

1 1 . 7803887784 

. 999990813211 

26 . 1 3 10601836 

12.8783237127 

. 999995507461 

27.8732276465 

13.866248646 

. 999997734746 

28.8 153850883 

15.0541735783 

. 999998825244 

3 1 . 6575625522 

16. 1420985125 

.999999374796 

Dbz  =  7. 32045232868 

Pdl  =  .5 

Pd2  =  .8 

Pd3  =  .99 

SIGMRL-TO-HOISE  PATIOS  CDB)  FOR  IHITIRL 

TRANSFORMATION 

7 , 18028480646 

16.5239686654 

26.8808804974 

12.6887488763 

21 . 2426978346 

31.503743108 

14.8466388857 

23.3091924052 

33.5584316626 

16. 268585645 

24.6986419313 

34.9428284825 

17. 3380310658 

25.7492651765 

35. 9905902914 

18.1 863823488 

26.5845353622 

36.8340131363 

18.812  0  6  7888  6 

27.3017783881 

37.5398641554 

SIGNAL-TO-NOISE  RATIOS  CDB)  FOR  NEW 

TRANSFORMATION; 

7.40601673703 

16. 6306689407 

26. 9721139815 

8.643746302  5  8 

18.3703967446 

28.6551266507 

10. 8741086655 

18.4515609092 

29.7168048237 

1  i  .  tt  t-  J  b  y  b  y  r'  'r* 

20.2655149832 

30.52121 18605 

12. 465801 1251 

20.93472161 1 1 

31. 1845531717 

1 3 . 0680388862 

21 . 5082203112 

31 . 7549390039 

13.5841676551 

22. 0147264365 

32.2573356581 

14.0624488581 

22. 4668270369 

32.7069765459 

14. 4847585558 

22 . 8760203338 

33. 1141652155 

14. 8684755601 

23.2498575243 

33. 4863332091 

15.2228081865 

23.5940083337 

33.8280709181 

15. 5485188303 

23.8128603223 

34. 1467104786 

15.8533503185 

24.2088871812 

34.4426853708 

16. 1373031587 

24 . 4878'^9623  3 

34.7197662984 

16.4  y  3  8  21  7386 

24.748160188 

34.9802228898 
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